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SUMMARY 

A theoretical study is made of the shielding 
mechanism whereby the rates of heat transfer are 
reduced near the stagnation point of two-dimensional 
and three-dimensional bodies when melting and 
vaporization occur simultaneously. The problem 
is approximated to its simplest form consistent with 
the mass-transfer and heat-transfer conditions . 
Simple results are presented which give the rate of 
mass loss at the body surface and the thickness of 
the liquid layer at the surface in terms of the rate 
of heat transfer to the unshielded body . 

INTRODUCTION 

The shielding of the nose of a body from 
excessive rates of heat transfer by melting or 
vaporization has been shown to be particularly 
attractive in view of the proportion of heat 
convected away from the nose in the boundary 
layer. (See refs. 1 to 4.) 

The extent of the shielding depends on the rate 
of mass loss, the manner in which the mass is 
removed, and its thermal capacity. Thus, it 
may he expected that the properties of the material 
in its molten and gaseous states, as well as in the 
solid state, play an important part in determining 
its suitability as an ablation shield. Since the 
rate of mass loss depends on the rate of heat 
transfer to the surface, the choice of material is 
very important inasmuch as this is the only 
control on the ablation rate left to the designer. 

Theoretical studies of surface melting have been 
made previously (refs. 1 and 2) where the convec- 
tion of heat by the molten material provides some 
shielding. When vaporization of the liquid is not 
considered, this shielding is necessarily limited 
to the heat convected by the liquid. When 


vaporization of the liquid is considered, however, 
the shielding can increase considerably because 
of the latent heat of vaporization and because 
the gaseous products now diffuse to the regions 
of higher (stagnation) temperature and arc con- 
vected in those regions (refs. 4 and 5). Whereas 
it is probably more desirable to have direct 
vaporization without a liquid film (since this 
precludes possible instabilities in the liquid film), 
it may be more efficient from overall heat-capacity 
considerations to have a material which, during 
ablation, produces both liquid and gaseous layers 
near the surface. 

This situation is considered in the present 
report by an extension of the simple' analyses 
made in references 4 and 5. An important 
feature of this more general type of ablation is 
the interaction of the gas and liquid layers due 
to vaporization at the gas-liquid interface. This 
interaction plays an important part in deter- 
mining the proportion of total ablated mass 
which is vaporized. Ablation with no vaporiza- 
tion of the liquid layer is obtained as a special 
case of this more general treatment. 

SYMBOLS 

x coordinate along solid surface 

y coordinate normal to solid surface 

with origin at solid surface 
y' coordinate normal to solid surface 

with origin at interface 
z transformed ^-coordinate 

Z arbitrary value of 2 outside gas 

boundary layer 

u component of velocity in z-direc- 

tion 

v component of velocity in ?y-direc- 

tion 
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Tii H-n 2 n x + n 2 


free-stream velocity in r-dircction 
constant in velocity distribution 
constant in velocity distribution 
for liquid 
temperature 

mass concentration of foreign gas 
effective mean concentration 
pressure 

pressure at stagnation point 
density 

coefficient of viscosity 
thermal conductivity 
coefficient- of binary diffusivity 
total rate of mass loss per unit, 
area of wall 

rate of vaporization per unit area 
of interface 

specific heat at constant pressure 
mean effective specific heat of gas 
mixture 

specific heat of liquid 
specific heat of solid 
latent heat of melting 
latent heat of vaporization 

fractions of gas and air in 1 mole of 
mixture 

molecular weight 
gas constant 

constant in Sutherland viscosity 
law 

, x dT 

JNusselt number, > r -j- 
T e — Ti Ay 

Reynolds number at interface, 
Schmidt number, — 

plJyi 

Prandtl number, 

skin-friction eo effi c i en t 
shear stress 

rate of heat transfer per unit area 
physical boundary -layer thickness 
t ra i \ sf ormed veloc i t y -houn darv- 
layer thickness 

trail sf ormed therm al -houn d ary- 
layer thickness 
transformed concentration- 
boundary-layer thickness 


A* 

D 

0 

Qo? Oh Q'i ? Qz 

J 



momentum thickness of gas 
boundary layer 

displacement thickness of gas 
boundary layer 
liquid-layer thickness 
t h ermal -1 ay or t h ick n ess with in 
solid 

enthalpy parameters (defined by 
eqs. (60a), (60d), (60c), and 

(60f), respectively) 
additional effect due to boiling 
(eq. 29(b)) 

functions of (defined by eqs. 
PiMf 

(00), (63), (64), and (68) or (69)) 
effective heat capacity (eq. (82)) 

external flow near stagnation point 

gas-liquid interface 

wall 

body 

liquid 

no ablation 

value at boiling point of liquid 

foreign gas 

air 

maximum 

two-dimensional 

axisymmetric 


THE SHIELDING MECHANISM 


The flow configuration considered is shown in 
figure 1. The coordinate system is fixed in the 
melting surface at the stagnation point so that the 
interior of the body moves with velocity the 
ablation rate, toward the surface. In the steady 
state there is a balance within the solid of conduc- 
tive and convective heat transfer (due to t lie 
motion of the solid). 

For sufficiently high rates of heat transfer a 
liquid film is produced by melting, the motion of 
which is controlled by the pressure gradient 
imposed by the outside stream and the shear 
stress at the interface. It is shown that, when 
no vaporization occurs, the presence of the liquid 
film has negligible effect on the air-boundary-Iaycr 
flow so that the interface shear stress assumes 
the value that holds in the presence of a fixed 
boundary. 

As the rate of heat transfer increases, the thick- 
ness of the liquid layer changes until the interface 
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Ficure 1. — Flow configuration. 


forces are in balance; the interface temperature 
increases to such a value that the convective 
shielding by the liquid layer is compatible with 
the rate of heat transfer at the surface required 
to produce the liquid. As the heat-transfer rate 
increases further, the rate of vaporization becomes 
sufficiently large to affect the character of the air 
boundary layer. This has two important effects 
on the interface conditions: Firstly, the rate of 
heat transfer across the interface is partly reduced 
because of shielding by the gas products and> 
secondly, the introduction of gas at the interface 
reduces the interface skin friction. As a result of 
this second effect the liquid layer tends to become 
thicker with an attendant reduction of gradients 
at the wall, although this trend is somewhat 


diminished and even reversed because the viscosity 
of the liquid decreases as the temperature increases. 
For a given rate of vaporization it is shown that 
the motion of the liquid film has a negligible effect 
on the gas boundary layer. 

In the gas boundary layer the gas produced by 
vaporization diffuses away from the interface and 
is convected as a mixture (of air and gas) parallel 
to the body surface. For exteme heating rates 
this convection can make the major contribution 
to the total shielding, since the gas is raised to the 
higher temperatures. This particular phase of 
shielding has been investigated in reference 4. 

Thus, it is seen that the ablation rate is con- 
trolled by the net heat transfer to the wall; this 
heat transfer, however, is itself greatly influenced 
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by the manner in which the products of ablation 
are removed from the wall. The analysis explains 
in a more precise fashion the balance of heat 
transfer and mass transfer which produces this 
shielding phenomenon. 

ANALYSIS 

HEAT TRANSFER IN THE SOLID 

For a solid moving with constant velocity 
toward the surface y — 0, the transfer of heat is 
governed approximately by the following equation: 

l( k *T\ 

dy \ dy ) 

Diffusion of heat 
to interior 


. dT 

mC ^ 

Convection of heat 
toward surface 


a) 


Equation (1) neglects temperature gradients 
parallel to the surface and is therefore a good 
approximation when the thickness of the solid 
ablat ion shield is large compared with the thickness 
of the thermal layer 0 within the solid (near the 
surface) in which the temperature changes rapidly 
from its value at the wall T w (the melting tempera- 
ture) to its value at the far interior TV (See 

%■ i) 

Equation (1) has the solution 


me by 

T^T b +{T u -T b )e~ (2) 

which satisfies the condition T- T w at ?/ — 0 and 
the condition T=T b as ?/— > — oo (with T w and T & 
constant). 

The thickness of the thermal layer near the sur- 
face is defined as follows: 


e= 



T—T b 
T ic — T b 



(3) 


In order to confine the high temperatures to a 
thin region near the surface, it is desirable to have 
a material of low conductivity k b and high specific 
heat c 6 . 

The rate of heat transfer away from the surface 
toward the interior at ?/= 0 is obtained from equa- 
tion (2) and is written 



=c b {T w —T b )m 


(4) 


where the subscript w— denotes that the surface is 
approached through negative values of y. The 
rate of mass loss m is unknown at this stage and 
must be determined by consideration of the liquid 
layer and the gas boundary layer. 

THE BOUNDARY CONDITIONS 

The method employed in the following sections 
assumes that the detailed variation of the velocity, 
temperature, and concentration profiles is of sec- 
ondary importance so that they may be approxi- 
mated in the simplest manner. It is important, 
however, that the boundary conditions which de- 
scribe the transfer of heat and mass at the melting 
surface and at the gas-liquid interface be carefully 
formulated and respected. 

Melting conditions .— 1 The heat-transfer condi- 
tion at the wall is 

(k ^ ) = L n m + c b (T w —T b )m (5) 

Heat transfer Latent heat Body heat 

to wall of melting 

The last term of equation (5) is given by equation 
(4) and is the heat required to raise the mass rh to 
the temperature T WJ in unit time, before melting. 

The mass-transfer condition is 


Pb^b — Pi v l,tc — M (fi) 

Solid lost Liquid produced 

from which the ablation rate v b is found when the 
rate of mass loss is known. 

A third condition is that which states that the 
liquid produced has no component of velocity along 
the wall, since at the wall the heat transfer, and 
therefore the melting, takes place along lines nor- 
mal to the wall; thus 

u t ,u>= 0 (7) 

Finally, 

T= T w = Constant (8) 

reflects the assumption that the melting tempera- 

ture T w is constant along the wall. 

Gas-liquid interface conditions. — The rate of 
heat transfer through the gas-liquid interface and 
the skin friction at the interface control the rate 
of melting at the wall and the manner in which the 
liquid is removed. These interface values, how- 
ever, are themselves influenced by the rate of mass 
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transfer across the interface and careful considera- 
tion of the boundary conditions is necessary. 

It is first shown that the motion of the liquid 
parallel to the wall has negligible effect on the gas 
boundary layer in general and on the interface 
values of skin friction and heat transfer in par- 
ticular. 

The pressure distribution imposed by the out- 
side stream is 


and, when the interface velocity is Ui %i , is 


TV 


-V=\ PeU >- 4 PCV 


(9) 


and the pressure in the liquid at the interface satis- 
fies the relation 

Pa— P^\ PlUl/^PlCfX* ( 10 ) 

Equation (10) is a statement that the decrease in 
pressure along the interface away from the stag- 
nation point exceeds the dynamic pressure ~PrW M 2 * 

because of the presence of viscous forces; the 
equality sign holds only when the viscous forces in 
the liquid layer are negligible. 

In order to match the variation with r, equa- 
tions (9) and (10) require that 1 


IpeC^lp'W 


(11) 


( 


dw\ U—Ui i 


(Here 5 is the thickness of the velocity boundary 
layer in the gas.) The difference in stress ex- 
pressed as a fraction of the value for no motion 
of the liquid is u lfi /U and therefore may be 
neglected. 

The effect on the rate of heat transfer across 
the surface is also easily estimated. The convec- 
tion of heat in the gas boundary layer parallel to 
the interface is given by 


pc p {T—T t )udy' 


If simple profiles 

T-T t y 

T'-Tt S 

u y' 

r~t 

are assumed when there is no motion of the liquid 
layer and an average value of pc p is assumed, then 
the convective heat transfer in the boundary layer 
parallel to the interface is 


bpCpiT.-T^U 


The velocity along the gas-liquid interface is, from 
equation (1 1 ), 

u lJ =C l x=-^U^(^) m U (12) 


When the interface velocity is u lti and the 
assumed velocity profile is 

. A 

U~U^V VJh 


For most liquids the density ratio p e ip t is less 
than 10“ 4 under the temperature conditions 
at which ablation can take place; therefore, 
nJV< 10" 2 . 

The effect of liquid motion on the interface skin 
friction is estimated very simply by the following 
approximations: The shear stress at the interface, 
for no motion of the liquid, is 



then the convective heat transfer is increased to 


(Hi 


1 \t 

The fractional increase is thus - -jj-7 that is, it may 

be neglected, and the corresponding reduction in 
heat transfer across the interface is negligible. 

The interface conditions can now be written 
without reference to the motion of the liquid. 
The mass-transfer condition is simply 


1 More precisely, it is the normal stress component that should he matched 

but at the interface this differs from p by a constant so that equation (11) is 

true exactly. 


PlVi.i — Pil'i — Wti (13) 

Liquid lost Gas introduced 
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The heat-transfer condition is 



Heat transfer Latent heat of Ileat transfer away 
to interface vaporization from interface 

A third transfer condition is that which states 
that no air crosses the interface; thus, 

-- < i5 > 

Diffusion of air Convection of air 
toward interface away from interface 

that is, there is a balance of diffusion of air toward 
the interface and convection of air, of concentra- 
tion 1 — W, away from the interface. 

The concentration of foreign gas also depends 
on the external pressure and is such that the 
partial pressure is compatible with the interface 
temperature. Thus if 1 mole of tin* binary 
mixture of air and gas at the interface consists 

0 f — Ul ... mole of gas and — ^ — mole of air, then 
nj+n 2 %+n 2 

the partial pressure of the gas is 
Pi_^ n x 

Pi ni + 7h 

and the mass concentration IE* is given by 

re 'it i 


The partial pressure p x and the interface temper- 
ature are related by the (simplified) Clausius- 
Clapeyron equation : 


a-AGHs) 

Pi 


where T B is the boiling temperature at the 
pressure p t . 

Elimination of Pi/Pt and rijn^ from the foregoing 
equations gives the following relation between 
Wt and T t : 


H',= 1 ('*' ( 7i r ^ _1 ) 


(1G) 


Iii addition, there are conditions for the con- 
tinuity of the components of stress at the inter- 
face; that for the normal component has already 
been given and yields 



( 17 ) 


Continuity of the tangential component of stress 
gives 

("!),-(«!)„ 

and, finally, there is no slip along the interface so 
that 

Ui=u iti (19) 


Equation (19) was tacitly assumed and used to 
show the negligible effect of liquid motion on the 
gas layer. 

When vaporization occurs, both m* and T t must 
be determined as part of the solution. If melting 
does not occur (that is, when the wall vaporizes 
directly with no liquid phase), then m — m*. 

The stress condition given by equation (18) is 
required to determine the magnitude of the liquid 
velocity u ipi which has an important effect on the 
liquid-layer shielding (even though it has a negli- 
gible effect on the gas layer). 

External stream conditions. — It is assumed that 
the following conditions hold outside the gas 
boundary layer: 

u u cA 


P=Pe 

ir -o 

These are the usual external conditions for stagna- 
tion-point flow; the last condition is a statement 
that there is no foreign gas outside the boundary 
layer. 


INTERFACE HEAT-TRANSFER AND SKIN-FRICTION 
RELATIONS 

Heat transfer. — Since the motion of the liquid 
layer has negligible effect on the gas boundary 
layer, the situation at the interface is exactly that 
which prevails when a solid surface vaporizes. 
The heat-transfer and mass-transfer results for 
such a problem arc given in references 4 and 5 and, 
therefore, only a brief summary is presented 
herein. 

The heat -transfer relations are found in terms of 
the Nusselt number N# Ut u Reynolds number R<, 
and Fraud tl number N Ptii which are defined at 
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the interface as follows: 


N„ u , i ~ 


x (<IT\ 

fXdy) t 


Z-I 


7? i= ^ J.2 

Mi 


Npr,i= : 


k] 


In addition, the Schmidt number is defined as 


The rate of heat transfer to the interface is then 
given by 

*-(* f T -> sb 

( 21 ) 

The following results (eqs. (22) to (29)) are 
taken directly from the analyses of references 4 
and 5 (where the subscript 0 denotes no vaporiza- 
tion) : 

For axisymmetric flow with constant pp, 

xt (tsH 0 ' 705 - 0 ' 065 0 - 0 } n -- ‘ 

(22a) 


and for two-dimensional flow with constant 


pm, 


-0.065 


o-m 


- 0.6 


(22b) 


( AV'Oo [ 0 ' 570 

For axisymmetric flow with variable pp, 

j vb (fbb 0 765 (j£T (M »> 


and for two-dimensional flow with variable pp, 

(§w) =0-570 (^Y'Vp,r 0 - 6 (23b) 


Equation (23a) is essentially that used in refer- 
ence 6 for an equilibrium gas with a Lewis number 
of unity. Results given by equation (23a) are 
compared with the exact results of reference 7 in 
figure 2, The effect of mass transfer by vaporiza- 
tion across the interface is given by 

504999 - 60 2 



Figure 2. — Variation of 


Abvu .i 

RX 2 


with 


p f p e 


PiUi 


Ax {symmetric 


stagnation point; no vaporization. 


I 1 \Nu.i\ / i 1 AT -0.g\ 

P,.i Ri in ~\N Pr . t Jt^'h V 3 A/ ‘ r -' ) 

[ i+ (b-0 Tr ] 


TYli 

(piHiC) u 


(24) 


where W is the mean concentration of foreign gas 
in the boundary layer and is given in reference 4 as 
a function of N Se ,i> Equation (24) is valid for 
both two-dimensional and axisymmetric flows. 
The last term in equation (24) represents the 
reduction in the rate of beat transfer due to the 
shielding effect of the mixture of air and gas in the 
boundary layer. A comparison of the results 
given by equation (24) with the exact results of 
reference 7 for the special case c Pi i=c p> 2 is shown 
in figure 3 for both variable and constant pp. 

The thickness of the velocity boundary layer 
with no vaporization S Ut0 is given in dimensionless 
form as 
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Figure 3. — Effect of vaporization on heat-transfer param- 
eter for two-dimensional flow. c Pi i — c v y, N Br>l — 0.7. 


where the superscript (1) denotes Np r%i — 1. The 
increase in bound ary -layer thickness due to mass 
transfer across the interface is given by 

oH — (2b) 

pjC 

Tlie thermal-bound ary-layer thickness 5 T and 
the concent ration-boundary-layer thickness <V flic 
given in terms of the velocity-boimdary-layer 
thickness as follows: 


N> 


0.6 . 


JYli 




sc ' f (p^O 


1/2 


/AWA ( \l 

\/^ l/2 /o ^ 


rrii 


( 28 ) 


3 (p,m,0 1/2 


From equation (28) it is seen that W t increases 
with t — ^ 7 wt 72 and reaches the value of unity when 

{PiPiL) 

CD 

{ ^ w a. * | 

Tfii 


/Vy a ,A (1 


(p^cy* v 0 . 6 i 


When TFj— 1, however, equation (16) shows that 
T t =T B) that is, the liquid at the interface has 
reached the boiling temperature. Since W t cannot 
increase above unity (by definition), equation (28) 
ceases to be true for 

/AWA a) 

rhi ^\lii l/2 Jo Wi.B 
(p^C)^ 2> v 0 . 6 1 (PiPiCy* 

^ Sc, i g 


The boundary condition given by equation (15) 


also shows that 


rrii 


/ and Wt arc not related 

(PiPiC) 112 * 

when Wi = l since equation (15) is satisfied for 
arbitrary values of m t . This is to be expected 
physically when the boiling temperature is reached. 

The heat-transfer relation given by equation 
(24) is modified when m <)> m n and is now ob- 
tained by the method of reference 5 as follows: 
The integral equation (13) of reference 5 is first 
written in the dimensionless form 


+(£-0£fI 


T.— T u 


T.-TtU 


Uz 


'N 


Pr , i 


n 


yu,i { ( C v, 1 1 \ 

, I/2 ) (PiPiC 


m i 


,cy /2 1 (PiUiCy 12 





The concentration of foreign gas at the interface 
I U is given by 


Substitution into this equation of II ',= 1 and 



N NuJ \ 1 m, 
n t 1 *) o + 3 (piUiCy/* 


(from eqs. (23), (24), and (26) of ref. 5) and the 
linear profiles given by equations (21) and (38) of 
reference 5 yields 
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1 N Nu j_ / 1 A x u ' 

N Pr . t 7?i' o 

where AT is defined by equation (40) of reference 5. 
Finally, substitution for K from equation (42) of 
reference 5 gives the form 

1 ^Vm,< / 1 A iV a A 

N PTti Ri'fi-VNpr.i lit 1 *) 0 

- J ( l _ 3 X "-- r " ") <p, r ,c>''’ <29a) 

where 

• 7 -G“- 1 I(- Vs - , 4) 


Tfli 


m t 

7W* 


0 1/2 

(29a) 


where 


(PiUiCy* /N Ku .i\n 

AL -°* 6 

m t VA7 /2 /o J 

^ Sc, i 


(29b) 

and represents the additional effect due to the 
boiling; the term is positive for c Pt \^>c Pt 2 and nega- 
tive for c Pi x<C.Cp t 2 . 

Skin friction, — The value of the skin friction 
whicli results when there is no mass transfer across 
the interface is first obtained by utilizing available 
exact solutions. 

Tim integral form of the momentum equation 
near the stagnation point is for axisymmetric flow 

and for two-dimensional flow 

which may be written 

t i ( jf ±\ /2 r j*t a ( u \ i 1 

sTWV L(PiVtO l/2 by \u)j(p^~C)^j 


where for axisymmetric flow 

A A* 

1= 3 

and for two-dimensional flow 


J= 2 -4-— 

In equations (31), 

W4 (-p)' 

A--fVi-SYh-fY: 


“-ro-5)— ro-?)- 

and 

«.= f f<v 

Jo Pi 

C yf p 

where z— — d?/' and Z is some arbitrary value 
Jo Pi ’ 

of ^ outside the boundary layer. The quantity 
A m is the momentum thickness and A* is the dis- 
placement thickness, based on a reference density 
Pi . The dimensionless ratios A m /$ w and A*/5 U , and 
hence /, are independent of the reference density, 
however. 

It is seen that I is a dimensionless quantity 
which depends on the details of the boundary- 
layer profiles; the method used herein to determine 
the effect of mass transfer on the skin friction 
(similar to that used in refs. 4 and 5 to determine 
the effect of mass transfer on heat transfer) avoids 
the evaluation of 7. Firstly, it is convenient to 
express equation (30) in terms of the Reynolds 
number and skin-friction coefficient 


••(I), 




2> iUi 


Equation (32a) can be expressed by use of the 
definition of Jit in the alternative form 


C f 7?J' 2 = 


Tj 1 

(p^CT' 2 Cx 


With these dimensionless coefficients equation (30) 
becomes 

(33) 
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The dimensionless thickness 8 U J already 

known and is given by equations (25) and (26). 

The skin-friction parameter Q C/7?j I/2 ^ is known 

from exact calculations (ref. 7) and was found to 
satisfy approximately the following equations: 
For axisymmetrie flow with constant pp, 


Q 6' / 7?, I/2 ) ( = 1.312— 0.489 (l -^) (34) 

and for two-dimensional flow with constant pp } 

Q <7,7? A 2 ) =1 .233-0.026 (l-O) (35) 

where the subscript 0 implies no mass transfer. 
The small variations with Prandtl number in this 
parameter have been neglected. 

A simple power law for the variation of the skin- 
friction parameter with pm was sought, similar to 
that for the heat-transfer coefficient (eqs. (23a) 
and (23b)) and valid in the range 


T e = 10,000° R 

CT? 

A 

A 

O 

7^= 2, 000° R 

(o.5<J<0.2) 


(The Sutherland viscosity law p °c 


2^3/2 

T+Su 


with 


Su— 200° R was used.) 

The following results were obtained from the 
exact solutions of reference 7 (and from inter- 
polation of these solutions): 

For axisymmetrie flow with variable pm, 


Q 5 (“)'■' (36) 

and for two-dimensional flow with variable pm, 

G w").- 102 (5ST (37) 


In figure 4 is presented a comparison between the 
results obtained by use of equations (36) and (37) 
and those of reference 7. 

The reduction of skin friction at the interface 
as a result of mass transfer across the interface 



I I L. . 1 L I 

0 .2 A .6 .8 1.0 

PePe 
PtH-i 

Figure 4. — Variation of skin-friction parameter with 

P*Pe tvt ■ *- 

No vaporization. 

PiPi 


is now found. It is assumed that, when mass 
transfer takes place, the momentum thickness 
and the displacement thickness arc increased 

in proportion to the increase in <5 U ; thus, and 

® u 

A* 

—> and lienee /, are assumed to remain unchanged. 

Ou 

Equations (26) and (33) then give 



CfEi 1 *) 

/o 


i. n 7? i/2 i 

2 t/Ui ^(p^cy' 2 


fr 


2 rhi 

‘(Pitney' 2 


(38) 


When expression (25) with the relevant equation 

/ p t \ 1/2 

of equations (22a) to (23b) is used for 5 Ut0 

and the relevant equation of equations (34) to 
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(37) is used for then equation (38) 

gives the following expressions for the skin- 
friction parameter \d f R^ 2 with mass transfer 
across the interface: 

For axisymmetric flow with constant pp, 

I C / 7?^ 2 = 1 . 312 ~0.489 (l-p) 


1.312-0.4891 


0.765-0.065 1 


and for two-dimensional flow with constant p/q 
\ C f Tt i ' /2 = 1 .233—0.626 (l-J-) 


m i 

(WtO ui 



1.233—0.626 ( 

-01 

0.570—0.065 ( 

- 0 ] 


(40) 


The results given by equation (40) are compared 
with the exact results of reference 7 in figure 5. 

The relations derived for constant pp (eqs. 
(39) and (40)) are presented only for the purpose 



Figure 5. — Plffect. of vaporization on skin-friction parameter for two-dimensional flow, pp constant. 
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of comparison with exact solutions as an ad- 
ditional check on the validity of the simplified 
approach used in this report. In the rest of the 
report only those expressions in which pp is 
variable are used. 

For axisymmetric flow with variable pp, 


and for two-dimensional flow with variable pp } 


I <7,7^/*= 1.02 



m t r 1 1.02 f Pe nx n , ~ 
(PiPiCy r2 l 3 0.570 Wi/ J w 


1 QR t w=i.2l 

2 



mt r 1 j 

(PiPiCy /2 L 3 


JL21 

0.765 



(41) 


A comparison of the results given by equation 
(42) with the exact results of reference 7 is shown 
in figure G. Again it is seen that the present 
method gives numerical results accurate enough 
for engineering calculations. 
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Now that the heat-transfer rate and skin fric- 
tion are known in terms of m t attention is turned 
to the liquid layer. 

FLOW IN THE LIQUID LAYER 

The mass balance in the liquid layer is written 
for the a xi symmetric case as 


velocity profile 


JL=1L 

Ui.t D 


(45) 


so that equation (43) for axisymmetric flow re- 
duces to 

m — iiii—CiPiD (4G) 


rh=2C,f Pi^-dy+mt 

J 0 

and for the two-dimensional case as 


(43) 


and equation 
reduces to 


(44) for two-dimensional flow 
m— 7» t =i CiPiD (47) 


m 



U 7 

pi —«y 

t 


Mass lost Mass convected 

by solid in liquid layer 


+ rhi (44) 

Mass transfer 
across interface 


The heat-balance equation in the liquid layer 
is treated in a similar way: 

For the axisymmetric case, 


In equations (43) and (44), D is the thickness of 
the liquid layer. 

The integral is evaluated by assuming a linear 


+2C, f Cipi{T—T K ) -^-dy+q a (48a) 
Jo U l> l 


and for the two-dimensional case, 


q, = L c rhi + 

Heat transfer Latent heat of 

from gas layer vaporization 


CiiTf—T^rhi 

Heating before 
vaporization 


+ C t ( D cMT-Tj^-dy + (48b) 

Jo Ui.i 

Convection of heat Ileat, transfer 

parallel to wall to wall 


where 

q w =m[L m +c b (T w — T b )] (49) 

The integral in equations (48) is evaluated by 
use of the assumed linear temperature profile 

T-J\ = ]L m 

TT-t w d 1 u; 

Insertion of equation (45), either equation (46) or 
(47), and equation (50) into equation (48a) or 
equation (48b) gives the following result : 

qz^Lv^i+CiiTi— TJrhi 

Tu) (m—Thi) + [L m -\-c b (T w — T b )]m (51) 

that is, the mass m—ihi removed as a liquid is 
raised through an effective temperature difference 
2 

of - (T t — T w ) during convection. Equation (51) 

o 

is valid for both two-dimensional and axis}mi- 
metric flows. 


The rate at which mass is convected away from 
the vicinity of the stagnation point in liquid form 
depends on the velocity at the interface u iti ~-CiX 
and the thickness of the liquid layer D as shown by 
equation (46) or equation (47). These quantities 
depend on the shear stress at the interface which 
is approximated by 



A comparison between the results obtained by 
use of the linear relation (52) and those obtained 
by using higher order profiles was made in refer- 
ence 1 and showed virtually no change in the 
shielding effect of the liquid layer. 

The relation between heat transfer across the 
liquid layer and the temperature difference is 
written in the approximate form as 

7 {T % -T w ) 1 , , , 

lit — -=2 (2m+SV) 


( 53 ) 
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which results from the assumption of a parabolic 
temperature profile. In equation (53) 

<li i — T w )mj-\-— cdTi 

— T w ) (m—mt)+qy, 

q w = [L m + c& ( T w T &)] th 

that is, q lti and q w are, respectively, the rates of 
heat transfer to the liquid layer at the interface 
and from the liquid layer at the wall. 

APPLICATION OF THE INTERFACE BOUNDARY CONDITIONS 

In the preceding sections expressions have been 
obtained for the shear stress and rate of heat 
transfer evaluated at both sides of the interface. 
The boundary conditions which relate correspond- 
ing quantities on either side of the interface are 
now applied. 

From equations (14), (21), and (54), 

c »^ t ~ T '''> (ww (a^~ iipO 

= IU +c l (T i -T„)]m i +*c t (T i - T m ) (m -m ( ) 



dimensionless parameters 


Nrr,i Ri l/2 


and i CjliV 2 


are already known in terms of through 

equations (23a) (or eq. (23b)), (24), and (41) (or 
oq. (42)). 

The rates of mass loss m and m t have been made 
dimensionless by the factor (piPiCy /2 in order to 
compare the present method with previous ex- 
act results. It is more convenient, however, to 
use the quantity (p e p e C) U2 for this purpose; as will 
be seen later, the use of ( p e p e Cy ,f 2 leads to some 
simplification of the analysis and presentation. 
Accordingly, equations (46) (or eq. (47)) and (55) to 

(58) are rewritten m terms of 

T u and D ( Pr ^) . Thus, from equations (55), 
(29), and (23a) (or eq. (23b)), 


where 


Qa ~ Ql {p2c)' n 

Qo— a iNp r ,r 0 ‘ Q Cp t2 (T c — T ( ) 


and, for the axisymmetric case, 


~h [Rm~\~Cb(T w — Tb)]th (55) 
and from equations (18), (32), and (52), 

(56) 

The concentration TT, is eliminated from equa- 
tions (16) and (28) to give 

Nnu A® ,1 rru 

Ri' /2 )o '3 (p.M.0 1/2 t ! 

A7 0 .6 Vli 

*'■* W,0 ,/2 

(57) 

when T„. 

Equations (53) and (54) are combined in the 
form 



D 2 ,Ui 




_ t _ [-^m4"Cii(r K — T b )]rii 


Thus, equation (46) (or eq. (47)) and equations 
(55) to (58) are five equations from which m, rii it 
T u D, and u, :i =C,x can be determined since the 


«,=0.765 (— ) 01 (60b) 

VPfMf/ 

and, for the two-dimensional case, 

„ 1= 0.57o(^) 01 (60c) 

VPiMi/ 

and 

<2, =c. b {T w - 7Y) +L m + 2 c, (T, - T m ) (60d) 

Q*=c t (T m - Tj+I^+dTt- T a ) (60e) 

Q a =L e + (c p +Jc p , 2 ) (1 ~\N Pr , r 06 ) {T t —T ( ) (60f) 

where 

Cp= z c P ,iW+ c Pi2 (l — W) (61) 

and, for T i <2T Bf 

J= 0 

It is seen that # 0 is associated with the rate of 
heat transfer to a nonablating body at a surface 
temperat ure T 7 *; Q x is the total beat absorbed when 

2 

unit mass of liquid is removed, -(T t — T w ) being 

u 
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the mean temperature rise, Q 2 is that absorbed by 
unit mass before vaporization, and Q z is that ab- 
sorbed during and after vaporization. 

Elimination of C t from equation (46) (or eq. 
(47)) and equation (56) with the aid of equation 
(41) (or eq. (42)) gives 


1— « 3 


m t 


<*2 


(PePtC) 


--K0T 


m — in i 

Tp^cW 2 


(62) 


where, for the axisymmctric case, 

« 2 = 1.21 (M*) 01 

\PiPi/ 

for t lie two-dimensional ease, 

a^O.ol^Y 1 

\PiPi / 

for the axisymmetric case, 


(63a) 


(63b) 




,1 _b21 (M A 0 

3 0.765 \PiPiJ 


1.21 


/pePe\° A 

\PiPiJ 


(64a) 


and, for the two-dimensional case, 

, 1 1-02 / peP ( \°' 


«3 = 


3 0.570 \PiHi / 


1.02 


9 / P'f* A 01 

~ \pipi/ 


(64b) 


Equation (58) is written by use of the definition 

A T Pl'iCl 

A'j — 


as 


1/2 ] 


TYli 


N 


Pr ^~Tc t 



Ipr.l.i L \A*l,i/ J V PeUe J 6 {p e H e C) ln 

n L m +Ci,{T w — 7\f] m _ . , 

and equation (57) becomes 


1 TO; /z»/i i\ \ 

1+ 3 =1+ m 2 ) 

AT 0.« VTj -^1 

5Cl< (P^.C)‘« 


(66) 


where is defined by cither equation (60b) or 
equation (60c). 


From equations (59), (62), (65), and (66) the 

four quantities J^yh’ T » “d 

1/2 

^ ) can be determined. 

Pi.i/ 

The detailed results obtained from these equa- 
tions depend on the numerical values of the 
thermal and diffusive parameters involved; how- 
ever, some general results are easily obtained 
upon examination of the equations. 

RESULTS AND DISCUSSION 

INTERFACE TEMPERATURE Ti 

It is noteworthy that for small values of 

- — l/2 r that is, for the lower range of T e , there is 

a single relation between the external temperature 
T e and the interface temperature T t which holds 
for both axisymmetric and two-dimensional flows. 

This is seen when - — , ... is made equal to zero in 

(p ef x e C) l/i 

equations (59), (62), and (65) and the quantities 


m 


(PePeCy' 2 

is 


ct 


and D ( — — ) are eliminated ; the result 

(PePe_y 
\PlPl,iJ 


c P . 2 (T<-T t ) 


Cb(T w — T &) - - Ci(Ti— T w ) 


AT 0.6 r~ 

x v Pt, i 

_Cb(T w 


Ci(Ti T w ) 


n 2 / 3 


TJ + Ln+y.iTi-TJ 


where, for the axisymmetric case, 


a 4 = 


= (1.21)^ /p.mA 0133 
0.765 \piPi ) 


and, for the two-dimensional case, 
(0.51) I/3 / n " \°- 133 


a 4 = 


:0.51)^ 3 (p e p A 01 

0.570 \PiPi ) 


(67) 


( 68 ) 


(69) 


The numerical factors in the expressions for a 4 
are both 1.4 (with an error of less than 1/2 percent) ; 

thus equation (67) with a 4 =1.4 ( ) is a 

\PiPi / 

single relation between T e and T 7 * since the quan- 


tities 


PePe 


and depend only on T t and T e 


PiPi,i PiPi 
Furthermore, in the range of high values of 
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T e — Ti where tlic condition L 

c b (T V} -T b )+L m +c l {T i -T w ) 


CvAT-Tt) 


applies ^but 


Tfii 


;rr ■ . is not neg 




; Y the terms 


. (fi^eCy* 

Q x and Q 2 may be neglected and equation (59) has 
the approximate' form, for Ti<^T B , 

e Ti) ^1 


Tv _/ \r o.e 

C p , 2 (TV— T iiNpr , r °* 6 V ^ *) 

- t ( N '- r - 0 (73) 

From equation (73) it is seen that it is not 
always possible for boiling to occur, for if 


\r o.e I 

— ^>SC,i rj 

3 

p ^ s 

L P, 2 AT 0,6 A 

^Pr,i q 


(74) 




m t 


{PeHeCY* 


(70) 


Equation (66) is written 
riii 


i+i 


3 a 1 (p l ./x f C) 1/2 . , Ml | ^ 

Uvf, 


(K) 


(71) 


y 0.6 m * 

axip'HjC)'* 


Tfii 


and elimination of the quantity — VwT/l g 1V09 

a relation between and Ti which is independent 
of a x and is therefore valid for both axisymmctric 
and two-dimensional flows; that is, 


Mo 


b(T_JA 

'ii \'i\ r B ) 


C v :i{T e Ti) l\ Prj 0-6 A f\ 


-1 INs, 


the right-hand side of equation (73) assumes nega- 
tive values and the equation cannot be satisfied; 
thus, equation (74) gives the condition under which 
boiling cannot occur -that is, when 

(a) Cp/Cp ' 2 is sufficiently large or 

(b) N Pr t is sufficiently large or 

(c) N Sc , t is sufficiently small. 

These conditions are precisely those which give 
maximum shielding by the 1 gas layer (ref. 4). 
It is noteworthy that the condition for boiling 
depends only on the thermal and diffusive prop- 
erties of the ablation material in the gaseous 
state; it is independent of the nose curvature, the 
pressure, the heat-transfer parameter, and the 
boiling temperature of the liquid. 

Figure 7 shows the range of N PTti and N Sc ,i 
in which boiling is possible. For any given value 


+( % - , -*-O-<7’.(- v - “-0 (72) 

(Even when the dependence of N Sc>i and N Pr j on 
Wi is considered, it may be shown that IF* itself 

771 r 

depends on - ^ M -C) 1/2 t ^ ir011 ^ 1 equation (28); 

thus, ax is eliminated with m { and equation (72) 
still does not contain cm) 

It is therefore expected that in general the inter- 
face temperature T\ (for a given material) depends 
only on the external temperature T e and is the 
same for both axisymmctric and two-dimensional 
bodies. 

LIQUID BOILING 

Another result of importance is obtained when 
it is noted that equation (72) gives the conditions 
under which possible boiling of the liquid occurs. 
As Ti^T B , equation (72) reduces to 


of boiling is possible in the region above the 

Cp j 2 

line but is not possible below the line. The effec- 
tive specific heat c p is given by equation (61) and 
the dependence of IF on N Prti and N Sc ,i is given 
in figure 8 of reference 5. In practice a desirable 
ablation material would have a high value of 
Cp/c Pt 2 and the boiling temperature would not be 
reached. 

When equation (74) is satisfied, equation (72) 
gives a limiting temperature above which 

the liquid temperature cannot be raised. Equating 
the right-hand side of equation (72) to zero (as 
T e —Ti becomes infinitely large) gives 



( 75 ) 


Dimensionless water-layer thickness, D 
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0 .5 1.0 1.5 2.C 

N Pr ,i 


Figure 7. — Conditions under which boiling is possible. 


This situation results from the shielding mecha- 
nism of the gas layer; each of the conditions (a), 
(b), and (c) tends to increase the shielding by the 
gas layer (as explained in detail in refs. 4 and 5) 
and, thus, restricts the amount of heat available 
for transfer across the gas-liquid interface. As 
the maximum temperature is approached and 
larger amounts of liquid are vaporized, any addi- 
tional heat made available by increasing T e is 
absorbed by the gas layer (since this layer must 
be heated in part to the high temperature T e ) 
which results only in additional convection of 
heat parallel to the gas-liquid interface without 
further increase in the heat transferred across the 
interface. 

MASS-FLOW RATES 

Asymptotic expressions for the rates of vaporiza- 
tion for large values of T e —Ti are also easily 



Fhtcjre 8.-— Variation of water-layer thickness with external temperature for axisymmetric flow. 
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obtained from equation (59) as 

riii 


a,.Vp r ,r 0 - 8 


(p e IX e C) 1 '' C_p / J \ N p -O-V 

c p . 2 \ 3 Pr •’ / 


(76) 


which depends only upon the gaseous properties 
of the ablation material. Thus, a relat ion of the 
following form is expected for large values of 
T-T t : 


r Wi 


L(p^<t) 1/2 J 

2— <?im ( a l)2 -dim 

r rii i 

1 

\_{p e p e cy 

J, 


0.745 (77) 


For small values of T e — T it where m t is neg- 
lected, equation (59) reduces to 


m 


a,Cp.2(7WMk£!^ 


(Pe» e cy /2 Cb(Ttc _ Tb)+Ln + 2 Cl (T-T m ) (78) 


and therefore, again, 


r m 

2 —dim (m) 2— dim 

Lw/vd 


r jh 1 

(«i)a 

A 

L (ftftO 

>/,>j 


0.745 (79) 


since there is a single relation between T e and T im 
The relations (77) and (79) can be expected to 
hold approximately in the entire range since they 
reflect the dependence of the mass-flow rates on 
the relevant heat-transfer coefficients. 

LIQUID-LAYER THICKNESS 

Since the relations (77) and (79) are approxi- 
mately valid in the entire temperature range, they 
lead to a similar ratio between the values of the 
liquid-layer thickness. Equations (65), (78), and 
(79) give the following result: 


Km 


(ai)2 


(«i)« 


K-)“1 

L \f*l r ts J 2 —dim 
As a check, equation (62) with 


=0.745 (80) 


rrii 


(p e M e C) 1/2 


=0 gives 


[" 

/ Cpi 
\Pi.< 

y/2- 

L__ 

[(£ 

y/2* 

], 

K 

c Pl v 

Mi ,*/ 

/2~ 

2 -dim t 

w 

12- 

2-dim 


=0.762 


The variation of the liquid-layer thickness 
depends largely on the variation of viscosity 
Pi * with temperature and on the rate of vaporiza- 
tion. For low heating rates, when 

m t ^ ^ m 
equation (62) shows that 



m 

mm"1 1/2 


Jp^eC ) 1 ' 2 

M 1 , 


For most liquids, however, the viscosity decreases 
as the temperature increases and this effect may be 
sufficient to cause D to decrease for part of the 
range of heat -transfer rates. When appreciable 
vaporization occurs, the interface shear stress 
imposed by the gas boundary layer is reduced 
(as evidenced by eqs. (41) and (42)); as a result 
the liquid layer adjusts itself by increasing the 

thickness until the interface stress y. ltt yy again 

balances that of the gas boundary layer. The 
detailed behavior of the liquid layer depends on the 
balance of these opposing effects — that is, the 
increase in the rat e of liquid flow and the reduction 
in shear stress which tend to increase the thickness 
and the reduction in liquid viscosity which tends 
to decrease the thickness. 

EFFECTIVE THERMAL CAPACITY 

A measure of the effective thermal capacity of 
the ablation material is given by 


where q ija is the rate of heat transfer to a non- 
ablating body at a surface tempera ture that is, 
from equations (21), (23a) (or (23b)), and (60b) 
(or (60c)) 


Qi , 0 — o: 1 (p e p £ ,(7) 1/2 c z , i 2(7 T e — Ti) Aiv, i °’ 8 

The parameter H eff is herein defined as the total 
amount of heat absorbed by unit mass of the 
material in solid, liquid, and gaseous states and 
should be as large as possible; 77 £// depends on 
T t — T t rather than on T e — T w although in practice 


STAGNATION-POINT SHIELDING BT MELTING AND VAPORIZATION 


19 


Equation (59) may thus he written 

fZco— \Cb{T w ~ T b ) + L m -\-~ ci{T t — T w )]m 

Ci(Ti — T w )-\ L v -\-(c p 

+Jc ,.*) (l ,Vp„r 0 ' 6 ) (T- (81) 

Equation (81) shows that in unit time an amount 
of heat q ii0 causes a total mass loss m of which m t 
is vaporized. The total mass loss m absorbs a 
quantity of heat c b {T w — T b ) + L m in the solid state 
2 

and - Ci(Ti~ T w ) as a liquid — that is, the average 
o 

2 

temperature rise in the liquid layer is -- (T t —T w )] 

o 

the amount w* is further raised through a tempera- 
ture difference i ( T t — T w ) to T t and absorbs an 

u 

amount of heat L v during vaporization and an 
additional amount (Cp+Jcp^) ^1— ^ A>,,r 0 ' 6 ^ (T c 

equations (59), (62), (65), and (66): 

T b =- 80° F 

7^=212° F 

(Pt = 1 atm) 

c 6 = 0.5 Btu/lb-°F 

144 Btu/lb 

A^pr, i — 0.7 

Mi- 18 


— Ti) as a gas. The latter contribution is the 
gas-layer shielding and because of its dependence 
on T e —T t is usually the most important effect. 

The effective heat capacity is then 

Heff—CbiT w T b ) Ci(T t — T w ) 

+[jj Ci(Ti—T w )-\-L v -\-(c p 
+Jc (i-J Vv.r°- 6 )(7;-r,.)] (82) 

Equation (82) shows that Tl cff increases rapidly 
with the rate of vaporization m-t since the dominant 
term is that which contains the factor T e —T t ; 
thus, a liquid of low boiling temperature is desira- 

ble in order to make -O as near unity as possible. 

APPLICATION TO A PARTICULAR EXAMPLE 

The results of the analysis are now applied to 
the melting and vaporization of ice; the following 
numerical values were used in the solution of 


Tjc—32 0 F 

7^ — 290° F 7 T * = 325° F 

(pt=4 atm) (Pi = 8 atm) 

c t — 1 Btu/lb-°F 

= 1070-0.5(7*, 

—32) Btu/lb 

A sc, f— 9.5 

M a =29 


When these numerical values* are used, TF is 
found to be 0.925 and "Cp/Cp , 2 is found to be ] .925; 
also, 

at 0.6 1 

Y=0.69 

tit n ft * 


Thus, the inequality of equation (74) is satisfied 
and boiling cannot occur; the limiting temperature 
at the gas-liquid interface is given by equation 


(75) as 

Ti imax —]9o° F (pi= 1 atm) 

T iimax = 265° F (p t = 4 atm) 

and 

T um IX =310° F (Pi=- S atm) 

Tlie relation between T e and T, was found to be 
the same for the axisymmetric and two-dimen- 
sional cases (with errors less than 1 percent in TV) 
for each of the values of p t . The ratios of dimen- 
sionless mass-loss rates and dimensionless liquid- 
layer thickness were approximately equal to 0.745 
as suggested by equations (77), (79), and (80); 



20 


TECHNICAL RETORT R— 1 0 — NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


these ratios were always in the range from 0.737 to 
0.757. In view of the relation between axisym- 
metrie and two-dimensional results only the 
axisymmctric results are shown (figs. 8 to 12). 

The variation of liquid -lay or thickness is par- 
ticularly interesting. As seen in figure 8, initially 

/ Pl C \ I/2 

the dimensionless water-layer thickness T) ( — — ) 


increases rapidly with external temperat ure T t and 
thereafter the thickness maintains an almost con- 
stant value, although the rate of liquid flow 
continues to increase (fig. 9). 

For values of p { less than 4 atmospheres the 
boiling point is low enough for the effects of 
vaporization and increasing mass flow to dominate 
that of the decreasing viscosity, and the water- 
layer thickness continues to increase slowly. At 
a value of p t of 4 atmospheres those effects just 
balance and the thickness remains constant at 
0.338 between 500° F and 14,000° F. At pres- 
sures greater than 4 atmospheres the boiling point 


is high enough for the decreasing viscosity to cause 
a local maximum before appreciable vaporization 
takes place; thereafter, the thickness decreases 
slightly and then slowly increases again as vapori- 
zation increases. For each of the three pressures 
considered the variation in water-layer thickness 
is small in the temperature range 1,000° F 
<7X14,000° F. 

The dimensionless rate of vaporization is vir- 
tually independent of the interface pressure p t as 
seen in figure 10. The total rate of mass loss, 
however, does vary with p L ) it decreases as p t 
increases (fig. 9). This effect is due to the increase 
in interface temperature with p t (fig. 11), which 
reduces the rate of heat transfer from the outside 
stream (proportional to T t — T { ) and also increases 

2 

the average temperature of the liquid ^ (T f — T w ) 

-| -T w . This increases the liquid shielding and 
reduces the heat available for melting. 

In figure 12 the effective heat capacity is seen 



Fig r re 0. Variation of total-rate-of-mass-Ioss parameter with external temperature at 
various values of p t for axisymmotric flow. 
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Fiufre 10. — Variation of rate-of-vaporization parameter with external temperature at 
various values of pi for axisymmetric flow. 


to increase with T e and also with the latter 
efTect is a result of the reduction in total mass loss. 
The greatest contribution to II e// is made by the 
vapor, especially in the higher range of external 
temperatures where the shielding by the vapor 
layer increases while that due to the water layer 
is virtually constant and makes a maximum con- 
2 

tribution ot -c t (T itM „—T w ) to TI eff . 


CONCLUDING REMARKS 

An approximate analysis has been made of 
steady melting and vaporization due to aerody- 
namic heating near the forward stagnation point 
of blunt bodies. 

The shielding mechanism of the liquid and gas 
layers which form over the body is discussed 
qualitatively and quantitatively and numerical 
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Figure II. Variation of interface temperature with external temperature at 
various values of p l for axisymmetric flow. 


results for a specific example, the melting and 
vaporization of ice, show that the greatest contri- 
bution to the shielding is made by the water-vapor 
layer when the heat-transfer rate is high enough. 

It is shown that, for any material which melts, 
the liquid layer cannot attain the boiling tempera- 
ture if the specific heat of the gas produced by 
vaporization lias a high enough specific beat- that 
is, if the condition 


Cp, 2 at o.b _ 

”3 

is satisfied (c p is the mean effective specific heat- of 
the gas mixture and the subscript i refers to the 
gas-liquid interface). In particular, when the 


Schmidt number N Sc is equal to the Prandtl num- 
ber A Tp rt this condition is satisfied if the specific 
heat of the gas produced c Pt t is greater than that of 
air c Pt2 . Thus, since Cp t C>c Pt2 is desirable in order 
to increase the gas-layer shielding, it is unlikely 
that the liquid layer will reach the boiling temper- 
ature when a desirable ablation material is used. 
A general conclusion is that shielding by the gas 
layer is much more effective than that of the liquid 
layer because of the large enthalpy difference 
across the gas layer. 

The rate of ablation is shown to depend on the 
behavior of the liquid and gas layers and upon the 
interaction of these layers through the conditions 
relating rate of heat transfer and shear stress at 
their common interface. In particular, the effect 
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External temperature, T e , °F 

Figure 12. — Variation of effective beat capacity II e // with external temperature T e for axisymmetric flow. 

of the rale of vaporization in reducing the inter- hand, tends to reduce the layer thickness. These 
face shear stress has an important effect on the effects are displayed in the results for the example 

behavior of the liquid layer. The rate of melting of melting ice. 

of the solid and the reduction in interface shear Lanqmt Rbsearch Center, 

stress cause the thickening of the liquid layer; National Aeronautics and Space Administration, 
the reduction in liquid viscosity, on the other Langley Field, Va., November IS, 1958 . 



24 


TECHNICAL REPORT R-10 NATIONAL AERONACTICS AND SPACE ADMINISTRATION 


REFERENCES 


1. Roberts, Leonard: On the Melting of a Semi-Infinite 

Body of Ice Placed in a Hot Stream of Air. Jour. 
Fluid Mech.j vol. 4, pt. 5, Sept. 1958, pp. 505-528. 

2. Sutton, George W.: The Hydrodynamics and Heat 

Conduction of a Melting Surface. Jour. Aero. Sci. 7 
vol. 25, no. 1, Jan. 1958, pp. 29-32, 36. 

3. Lees, Lester: Convective Heat Transfer With Mass 

Addition and Chemical Reactions. Presented at 
Third Combustion and Propulsion Colloquium of 
AGARD (Palermo, Sicily), Mar. 17-21, 1958. 

4. Roberts, Leonard: A Theoretical Study of Stagnation- 


Point Ablation. NASA TR R-9, 1959. (Supersedes 
NACA TX 4392.) 

5. Roberts, Leonard: Mass Transfer Cooling Near the 

Stagnation Point. NASA TR R-8, 1959. (Super- 
sedes NACA TN 4391.) 

6. Fay, J. A., and Riddell, F. R.: Theory of Stagnation 

Point Heat Transfer in Dissociated Air. Jour. 
Aero. Sci., vol. 25, no. 2, Feb. 1958, pp. 73 85, 121. 

7. Beckwith, Ivan E. : Similar Solutions for the Compres- 

sible Boundary Layer on a Yawed Cylinder With 
Transpiration Cooling. NACA TN 4345, 1958. 


U.S, GOVERNMENT PRINTING OFFICE; I960 


